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On the Numerical Solution of Linear Integral Equations. 

By H. Bateman. 

(Communicated by Prof. E. T, Whittaker, F.E.S. Received October 29, 1921.) 

§ 1. Introduction. 

The numerical solution of linear integral equations of the types studied by 
Volterra has formed the subject of a recent memoir by E. T. Whittaker.* 
Numerical methods are needed also for the solution of the linear integral 
equation studied by Fredholmf and HilbertJ and the object of this paper is 
to describe a method which may sometimes be useful. 

The linear integral equation of the second kind, 

f(s) = (/> (s)_\ j\ (s, t) tj> (0 dt, (1) 

may be solved for the unknown function </>(£) either directly or indirectly. 
In the direct methods of solution the function $(t) is expressed by means of 
infinite series of terms involving repeated integrals. In the so-called method 
of Neumann only one infinite series is used, while in the more complete 
method of Fredholm the function cj)(t) is expressed as the ratio of two infinite 
series which converge for all values of the parameter X. The method of 
Neumann has been employed on many occasions to obtain numerical results 
and is generally more convenient to use than Fredholm's method on account of 
the great complexity of the expressions occurring in Fredholm's series. There 
are occasions, however, when the Neumann series fails to converge, or converges 
only slowly, and then some other method such as Fredholm's must be used. 
It seems desirable if possible to devise simple approximate methods which 
possess some of the advantages of Fredholm's method, because in many cases 
the evaluation of repeated integrals becomes very tedious and the use of the 
direct methods of solution becomes impracticable except for a rough approxi- 
mation. 

Of the indirect methods of solution, we shall be concerned here with one 
which depends on the use of an approximate expression for the kernel of the 
integral equation. E. Schmidt§ has already sketched a method in which the 
kernel k(s, t) is replaced by an approximate expression 



n 



tc (s, I) = t a m (s) /3 m (t), 

in = 1 

* E. T. Whittaker, ' Roy. Soc. Proa,' A, vol. 94, p. 367 (1918). 

t I. Fredholm, ' Acta Math.,' vol. 27 (1903). 

| D. Hilbert, « Gott. Nachr.,' 1904. 

§ E. Schmidt, ' Math. Ann./ vol. 63, p. 433 (1907) ; vol. 64, p. 161 (1907). 
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consisting of a finite number, n 9 of products of functions of s and t. It is 
known that the solving function for this approximate kernel can be written 
clown in finite terms* and the solution of the integral equation with this 
approximate kernel can be regarded as an approximate solution of the original 
integral equation. The success of the method evidently depends upon an 
appropriate choice of the functions a m (s), /3 m (t). Schmidt has something to 
say on this matter; in particular, he indicates how the functions u m (s), /3 m (t) 
must be chosen so as to make the double integral 



o Jo 



[fc (s, t)—7c (s, t)Y ds dl 



urn. The required functions a m (s) t f3 m (t) turn out to be the first n 

normal orthogonal functions associated with the kernel tc(s, t) in a 

manner described by Schmidt. These functions are generally difficult to find, 

consequently this best method of approximating to the kernel is generally 

useless for an easy numerical solution of the integral equation. 

In the following paper we shall show that by making a different choice of 
the functions u m (s) t j3 m (t) we can obtain formulae which are quite suitable for 
numerical work. 

§ 2. Method of Approximation. 

We shall use the term solving function to denote a function K(a, t) such that 
the solution of the integral equation (1) may be expressed in the form 

0(s)=/(s) + \fK(M)/(O* ( 2 ) 



when f(s) is any continuous function. Our method is now based on the 
following general theorem.f 

Let k(s, t) be a kernel for which the solving function K(s,t) is known, then 
if/i(s), / 2 (s) ...f n (s)'}'gi(t)> 9$) ...^»(0 are given continuous functions and 
aim (I, m = 1, 2, . . . n) given constants, the solving function H (s, t) for a 
kernel h (s, t), defined by the equation 



k(s, t)—h(s t t) /i(s) 
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* E. Goursat, 'Bull. Soc. Math. France,' vol. 35, p. 163 (1907) ; H. Lebesgue, ibid., 
vol. 36, p. 3 (1908) ; E. Schmidt, loc. ciL 
+ < Mess, of Math., 5 April (1908). 
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is given by the equation 

K(M)~H(M) ft 00 <h(*) 
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%m (0 = 0« (0 + M #™ ( S ) K ( S > *> > • 

Aj m = «*,» + X I # («) <£ m (a?) da? 



(5) 



The functions occurring in this theorem are all supposed to be completely 
defined in the intervals (0=s^l),(Q^t^l) and to be continuous throughout 
these intervals. 

We are interested chiefly in the case when tc (s, t) = 0, K (s, t) = 0, then 

<M 5 ) == ,/K s )> %m(0 = $w(0> an d the value of A im is at once known. 
Let us suppose that w(s, t) is the kernel of an integral equation 



*(s) = /3(s)~\ [ w(s, t)fi(t)dt 



(6) 



whose approximate solution is desired. 
Putting 

fm (s) = —wfa i m ), gi(t) — —w(s h t), a im = — w(s z , * m ), (7) 

where s h s 2y ... s n ; £ b £ 2 , ^ are particular values of 5 and t belonging to the 
interval (0, 1) it is clear that equation (3) becomes 

h (s, t) w (s, t{) vj (s,t 2 ) to (s, t n ) I = 0. (8) 

W (Si, W (Si, ti) w (s h t 2 ) w (s h t n ) 

W (S 2 , t) '10 (S 2} ti) W ($ 2 , t 2 ) W (S 2 , t n ) 



W (S n , t) W (S n , t X ) W (S n , t 2 ) V) (s n> t n ) 

This equation defines a function h (s } t), which satisfies the conditions 

A (s, t m ) = w (s, t m ) (0 ^ s ^ 1) (m = 1, 2 . .. n) 

h (s h t) = w (s h t) (O^sl) (J = 1, 2 ... ra) 

The function h (s, t) thus agrees in value with the function w (s, £) over a 
network of lines in the (s, plane ; it thus represents a good approximation 



(9) 
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to the function w (s, t). It follows, moreover, from formula (4) that the 
solving function H (s, t) can be written down, and so an approximate solution 
of the integral equation (6) is furnished by the formula 

j3 (s) = a (s) + X [ H (s, t) a (t) dt. (10) 

It should be noticed that the above approximate representation of w (s, t) 
becomes exact when w (s, t) is of the form 

■iv (s> t) = % p m (s) q m (t), 

m = 1 

it fails when w (s, t) is the sum of less than n products of type p (s) q (t), for 
then the co-factor of h($, t) in the determinant (3) vanishes. Approximate 
values of the singular values of X for the integral equation may be obtained 
by equating to zero the coefficient of H (s, t) in equation (4). It should be 
remarked that, when n is even, the approximate values of X thus obtained 
may be complex even when the exact values of X are all real. 

§ 3. Illustrative Example. 

Let us apply our method to Hilbert's integral equation, in which the 
kernel w (s, t) is defined by the equations 

s(l — l) s^t 
w(s, t) = 

I (l-~s) . s^t 

We have 

f 1 
Mm = — w (s h t m ) + X iv (s h x) w (x, t m ) dx. 

Jo 

Now w (s. x) w (x. t)dx = 

Jo ~yt(l-s)(2s-s 2 -~t 2 ) s^t, 

hence the value of Ai m is easily found. 

We shall interest ourselves simply in the way in which our method 
enables us to approximate to the singular values of X, for which the 
homogeneous equation 

(f) (5) = X iv (5, t) <fi (t) dt 

Jo 

can be satisfied. It is known that the exact singular values of X are 
7r 2 , (2tt) 2 , (37t 2 ), ..., while, according to our method of approximation, they 
should be given approximately by the determinant | Aj m | = 0, the degree of 
approximation being better the larger the value of n. Let us take n = 5, 

si = h = £, s 2 = t 2 = h s s = h = h ^ = ^ = i, So = ^5 = & 
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then if X = 216 a? the equation for X may be written in the form 
5(l-lte) 4(1 -19a;) 3(1 — 26a?) 2(1 — 31a?) 1-34® 
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3(1- 
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4(1- 


-19a?) 
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-10a?) 



= 0. 



The determinant is easily evaluated by repeatedly subtracting rows and 
columns and we obtain the equation 

130a? 5 -441a?* + 488a? 3 -206a? 2 + 30a?-l = 0, 

■i.e., (a?-l)(2a?-l)(5a?-l)(13a? 2 -22a? + l) = 0. 

Approximate singular values of X are thus 

\ x = 10-09, X 2 = 43-2, X 3 = 108, X 4 = 216, X 5 = 355*2. 

The first is equal to tt 2 within about 2 per cent,, the second differs from 
4ir 2 by about 9 per cent, while the errors in the other cases are larger still 
If with n — 3 we take 

Si = ^ = | r S 2 = h = i, 5 3 = h = f , 
we shall find the approximate values 

Xx = 10-37, X 2 = 48-13, X 3 = 12613. 
With n = 1 and s x = ti = J we find simply the approximate value 

In the case when w = 5, the percentage error in the value of X is roughly 
proportional to X. When we are dealing with the solution of the equation (6) 
instead of the homogeneous equation, we can generally expect the percentage 
•error to increase with X. If we wish our answer to be correct within 1 per 
cent when the value of X is in the neighbourhood of the first singular value, 
we should take n = 7 at least, with a value of X equal to about half the first 
singular value, it may be sufficient to take % = 5, for values of X larger than 
the first singular value, still higher values of % such as 9 or 11, may be 
.needed to give the desired degree of approximation. 

§ 4. Polynomial Kernel. 

The method of § 2 can be extended without difficulty to the case of a 
system of integral equations by making use of Fredholm's well known device. 
It may also be extended in other directions ; we may, for instance, write 

fm 0) = fC (S, t m ) — W (S, t m \ gi (t) = K (S h t)—W ($ h t), 

ai m = K (S h t m ) - W (S h t m ), 
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where tc (s, t) is a kernel for which the solving function K (s, t) is kno wn.. 
The function h(s, t) is then equal to iv(s, t) — k{s 3 i) over the lines of a 
network. 

Another method depends on the use of an approximate polynomial 
representation of the kernel w (s, t). Let us now put 



then 



k (s, i) - 0, f m (s) = s» g m (0 = ** K (s, t) 

A im =: a lm "4" ^ 
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If now the expression for h (s, t) is 

h (s, t) = — 



■5'i 



/J 



= m = 



where & M r = 0, m + r > % 

we may determine the constants a hn by solving tfre equations 
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It is clear from these equations that a mr = 0, m + r < w-fl, also 
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Our object now is to obtain an approximate representation of a kernel 
ye (5, t) by means of a polynomial in s and £, or by some other double series of 
functions of s and t. 

Some interpolation formulae for functions of two variables have been 
studied by S. Namuri * If we write 

A x l f(x + h, y) — A x l f(x, y) = Aj +l f(x } y) 9 

A^Ajf(x, y + k)-~A y ™A z l f(% } y) = A y » +1 A^/fo 3/), 

* S. Nanrnri, ' The Tolioku Mathematical Journal, 5 vol. 18, p. 309 (1920). 
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for all values of 1, m 0, i, 2, ... , it is easy to see that the operators A^ and A y . 
are permutable, in other words we have the relation 

A/ 1 A**/(*, ?/) = &J & y m f{x, y) 

for all values of y. Introducing the notation 

's\ _ s(s-~ 1) ... (5— fc+1) 

,&/ 1 7 2 .771 ' 

the desired interpolation formula may be written 

f(a + sh, b + tk) = "s 2 1 (j) (' ) A^ A/7>, S) + R M „ 

x = a-\-sh, y — b + tlc, 
where the remainder K^ has the form 

;) »&;/i*>»)+{S) * $f(*. «-(;) (t) ^ g^/(-. *V 

where « lies between the least and greatest of the quantities 

a, a-\-h, a-\-2h ... , a + (/x — l)h } x, 
and ft lies between the least and greatest of the quantities 

b, b + k, b + 2k ... , & + (v — l)&, 2/. 

If the function /(#, 2/) is actually a polynomial of degree /x — 1 in x and of 
degree v — 1 in y, the remainder vanishes. 

The above method of interpolation is suitable for the case in which the 
function f(x, y) is known at the points of intersection of two sets of lines 
parallel to the axes of co-ordinates. Mamuri has extended the formula for 
the case in which the interval between two successive lines is not constant, 
but it is desirable to have a formula which is applicable when the points at 
which f(x, y) is known are distributed at random. 

Let us consider the simple case in which there are six points, Pi, P 2 , P 3 , P 4 . 
P 5 , P 6 , and let G r = be the equation of the conic through all the points 
except P r , also let c r be the value of G r when the co-ordinates of P r are 
substituted in C r , then if p r denotes the value of f(x, y) at the point T r , the 
required formula is 

/ (^ y) — 2 Pr -7. 

r = 1 ^r 

When the six points are the six intersections of four lines Li = 0, L 2 = 0,, 
L 3 == 0, L 4 = 0, we may write 

Ci = L2L3, G2 = L3L1, C3 = LiL 2 , C4 = LiL 4) C 5 = L2L4, Gq = L3L4. 

These two methods, which are well known, may be extended, for when 
there are %n(n+l) points at which the value oif(x,y) is known, we may use 
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a curve of the (?i~-l)th degree instead of a conic, and when the %n(n-\-l) 
points lie at the intersection of n + 1 lines we may represent each of these 
•curves as the product of (n — 1) lines. 

§ 5. Possible Applications, 

It is known that many problems in the theory of elasticity may be reduced 
to an integral equation of the second kind or to a system of integral equations 
of the second kind.* In particular, Fredholm has indicated a method of 
reducing a two-dimensional elastic problem to a pair of integral equations 
with kernels that do not become infinite within the range of integration. 
The torsional problem in the theory of elasticity! and many electrical 
problems are equivalent to well known potential problems that are reducible 
to integral equations of the second kind. Certain problems relating to the 
flexure of beams are likewise reducible to integral equations. 

It is thought that the method of § 2 may be sometimes useful in dealing 
with such problems. It is true that the work will be laborious, but this 
seems. unavoidable when a numerical solution is desired. 

Problems relating to the critical strength of struts and structures, the 
vibrations of rods and membranes, the whirling of shafts, etc., may also be 
treated by this method, the differential equations being replaced by integral 
equations. It is thought, however, that in these cases the methods of approxi- 
mation proposed recently by other writersj may be more convenient to use. 

In the case of the equation of a vibrating membrane 

d-xr oy 
we may employ Whittaker's solution§ 

V = ' e ik ( xcma+ ^ sin ^<p(a)cloc. 

Jo 

* See, for instance, E. and F. Cosserat, ' Comptes Rendus,' Paris, vol. 133 (1901); 
I. Fredhohn, ' Arkiv for Math, och Fys.,' vol. 2 (1906) ; G. Lauricella, ' Atti Lincei ' (5), 
vol. 15 l5 pp. 426, 610, vol. 15 25 pp. 75-83 (1906) ; ' II Nuovo Cimento ' (5), vol. 13 (1907), 
four papers ; ' Rend. Lineei ' (1907) ; ' Acta Math.' (1908-9) ; T. Boggio, ' Atti Lincei, 5 
vol. 16 2 , pp. 248, 441 (1907); Mareolongo, 'Toulouse Ann.' (1908); J. Hadamard, 
'Memoires de l'lnstitut ' (1908); A. Korn, 'Ann. Ecole Normals ; (3), vol. 25(1908); 
Zaremba, ibid., p. 337; A. Haar, 'Gott. Nach. ' (1907); O. Tedone, 'Rend. Lincei,' 
p. 232 (1904) ; C. E. Weatherburn, 'Phil. Mag.,' vol. 32, p. 15 (1916) ; ' Mess, of Math./ 
vol. 47, p. 114(1917). 

t Of. L. Bairstow, ' Roy. Soc. Proc.,' A, vol. 95, p. 457 (1919). 

X W. L. Cowley and H. Levy, « Roy. Soc. Proc.,' A, vol. 95, p. 440 (1919) ; 'Phil. Mag.,' 
vol. 41, p. 584 (1921); A. A. Griffith, 'National Advisory Committee for Aeronautics,' 
R. and M. 545 (1918-19) ; R. V. Southwell, 'Phil. Mag.,' vol. 41, p. 419 (1921). 

§ 'Math. Ann.' (1.903); "Whittaker-Watson, "Modern Analysis." 
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If the boundary condition is 

V = when x = r(0) cos 0, y = r(0) sin 0, 

we obtain the integral equation 



<aTT 



J 

which may be regarded as a limiting case of the homogeneous equation 

f27T 



We may apply our method of approximation to obtain the possible values 
of the functions \x (k). Equating these to zero we obtain the possible values 
of h 



The Atomic Process in Ferromagnetic Induction. 
By Sir J. Alfred Ewing, K.C.B., F.RS. 

(Beceived November 29, 1921.) 

The object of this paper is to amend, in an important particular, the 
theory of ferromagnetic induction put forward by me more than 30 years 
ago,* and to describe a new model. That theory was itself a modification 
of the earlier theory of Weber. To Weber is due the fundamental notion 
that a substance contains minute particles, each of which acts as a magnet, 
and that in the process of magnetising a ferromagnetic substance these 
are turned into more or less complete alignment. The ultimate magnetic 
particles used to be called " molecular magnets " : we now recognise them 
as attributes of the atom, not of the molecule, and (in all probability) they 
derive their magnetic moment from the circulation of electricity in electron 
orbits or in ring electrons. What turns is not the molecule nor the atom, 
but something within the atom. 

The characteristics which distinguish ferromagnetic substances from other 
paramagnetics are: (1) the much larger amount of magnetism they can 
acquire under the action of an impressed field ; (2) the fact that the acquired 
magnetism tends towards a saturation limit when the field is progressively 
increased ; (3) the fact that the acquired magnetism shows hysteresis with 

* " Contributions to the Molecular Theory of Induced Magnetism," ' Roy. Soc. Proc, ' 
vol. 48, p. 342 (1890) ; * Phil. Mag.,' September, 1890. See also ' Magnetic Induction in 
Iron and other Metals/ Chap. XI. 
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